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ABSTRACT 

A hydrodynamical description for vortex states in type II superconductors is presented 
based on the time-dependent Ginzburg-Landau equation (TDGL). In contrast to the fa- 
miliar extension of a single vortex dynamics based on the force balance, our description 
is consistent with the known hydrodynamics of a rotating neutral superfluid and correctly 
includes informations on the Goldstone mode. Further it enables one to examine non- 
local conductivities perpendicular to the field in terms of Kubo formula. Typically, the 
nonlocal conductivities deviate from the usual vortex flow expressions, as the nonlocal- 
ity parallel to the field becomes weaker than the perpendicular one measuring a degree 
of positional correlations, and, for instance, the superconducting contribution of dc Hall 
conductivity nonlocal only in directions perpendicular to the field becomes vanishingly 
small in the situations with large shear viscosity, leading to an experimentally measurable 
relation p xy ~ p xx among the resistivity components. Other situations are also discussed 
on the basis of the resulting expressions. 
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§1. Introduction 

At present, understanding theoretically the nonlocality 1 ' 2 ) of linear resistivity in a 
clean system (with no pinning) seems to be one of central problems on the vortex states of 
a type II superconductor. By analogy to the usual viscous fluid, Marchetti and Nelson 1 ) 
have previously argued within a linear hydrodynamics that, deep in the liquid regime, a 
shear-viscous force proportional to — dj_dts will take the place of the shear term in the 
elastic force f e i in the following force-balance equation for the vortex lattice 

9 t (-r lS + r 2 (zxs)) + JxB = f el) (i) 

and that consequently, a positive nonlocal term proportional to q\ should exist in the 
conductivity perpendicular to B. In eq. (1) (and through this paper,) s denotes the dis- 
placement field of vortex positions, J the external current, B the uniform flux density 
parallel to z axis, q± the wavevector perpendicular to B, and i\ and T 2 will be given later. 
This equation is essentially an extension 3 ) of the single vortex dynamics to interacting 
vortex states. Subsequently, the presence of other nonlocal terms was proposed 2 ) neglect- 
ing any sample disorder, although experiments have been performed in heavily twinned 
(disordered) samples 4 ), and, the presence deep in the liquid regime of a large conductiv- 
ity nonlocal in the field direction was argued. It is important to theoretically clarify to 
what degree these proposals can be justified beyond the phenomenology 1 ' 2 ' 5 -* and within a 
basic dynamical model for type II superconductors such as the time-dependent Ginzburg- 
Landau equation (TDGL) which correctly describes the uniform linear dissipation in terms 
of Kubo formula. 6 ^ Actually, the eq.(l) is not correct in the following senses. Firstly, it 
is not compatible with calculations of conductivities based on the use of Kubo formula. 
Secondly, it is not clear from eq. (1), in which f e i is independent of J and of the time 
derivative dt, how the possible viscous (nonlocal) terms in the dc conductivity are changed 
as the vortex lattice is formed. Further, the dispersion ~ g 4 of the shear mode found as 
a phase fluctuation formally and generally 7 ' 8 ) cannot be seen in eq.(l). Improving this 
is important in order to avoid misunderstanding conclusions on the phase correlations in 
refs. 7 and 8. 

In the present paper, a correct hydrodynamical approach, or equivalently, a dynamical 
harmonic analysis on the transport in vortex states, particularly in the vortex lattice, 
is presented on the basis of TDGL. Our approach has no difficulties mentioned above 
and makes it possible to examine nonlocal conductivities above and below the melting 
transition, and several consequences of the resulting conductivities are discussed. As has 
been understood through studies 9 ) on nonlocal and static linear responses, the use of 
the harmonic analysis for 3D and layered systems often leads to misunderstandings on 
physics in the liquid regime, and hence, we will not consider situations with negligibly 
small positional correlations of vortices. 

§2. TDGL hydrodynamics 

Our starting point is TDGL describing dynamics of the order parameter ip in an 
isotropic type II superconductor 

0(7 + iVW + [&(|Vf - Po) + <o(-i<9 " = °; ( 2 ) 
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where curlA = Bz, 7, a, and b are positive, 0o is the flux quantum, £ the coherence 
length, and the po is the spatial average of the mean squared order parameter. The inclu- 
sion of sample anisotropy is straightforward and will be done later. Since the linearized 
version of eq. (2) is considered through the paper, the noise term introducing fluctu- 
ations in equilibrium does not have to be included, and growths of time scales will be 
phenomenologically included later as the only fluctuation effects. Within the harmonic 
analysis given below, spatially varying gauge fields other than the external disturbance 
(see below) have to be excluded in order, for the derivation of conductivities in vortex 
states, to become consistent with the corresponding one, for instance, of the ac conductiv- 
ity in the Meissner state Cij(u>) = irp s oS(u)Sij valid irrespective of the details of dynamics, 
where p s o = 2(27r£o/</>o) 2 Po° is the mean field superfluid density. 

Particularly upon calculations of linear responses, it is convenient to work rather 
within the harmonic approximation of the corresponding Euclidean action with gauge 
field disturbance S A (and with h = 1) 
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Since we are interested only in frequencies and wavevectors accompanying 8A and sum- 
mations with respect to them are not performed, calculations of conductivities based on 
the action (3) are equivalent to the linearized analysis on eq. (2) and hence, formally 
independent of temperature. Here r is the imaginary time, is the temporal Fourier 
transform of ip(r), O is Matsubara frequency, (3 the inverse temperature, and ipo the mean 
field solution of ip. The superconducting (and real) part o~ij of dc linear conductivity tensor 
should be always calculated in terms of Kubo formula 
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where fio is the external (Matsubara) frequency, i, j = x or y, and the absence of dc 
uniform (q = 0) superfluid density was used (see below). More generally, the ac nonlocal 
conductivities are defined as the imaginary parts of the response function p S)ij (q, cu + i0 + ) 
divided by the real frequency 10. However, since general expressions of the ac nonlocal 
conductivities are complicated and they always become the vortex flow results in q — > 
limit ( of the present formulation ), the dc nonlocal conductivities will be only considered 
below. Since the harmonic modes appearing in eq. (4) are spatial and temporal variations 
around the mean field solution linearly excited by the gauge field disturbance, the linear 
response quantities in vortex states resulting from such a harmonic analysis 10 - 1 are not 
accompanied by the temperature In this sense, this approach is not a fluctuation 

theory and is invalid for the linear dissipation parallel to B which is a consequence of 
critical fluctuations. 6 ) On the other hand, the vortex flow conductivities are essentially 
mean field results. Therefore, we focus below on the configuration perpendicular to B. 
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Before presenting general results, it is useful to first see results on the linear conduc- 
tivities in the Landau level (LL) approach 7 ) within the lowest and next-lowest LLs. In such 
a high field approximation the mean field solution and the shear elastic mode are found 
within the lowest LL, and the next-lowest LL, which does not participate in constructing 
the mean field solution, 7 ' 9 ) provides the uniform displacement of vortices (see eqs.(4.12) 
and (4.16) in ref.7). Since, in constructing basis functions of eigenmodes, Matsubara fre- 
quencies in the Euclidean action play similar roles to the wavevectors q z parallel to B, 
including dynamical terms is easily performed following ref.7 . For later convenience, an 
amplitude-dominated mode in the lowest LL, denoted by dp, will also be included. Then, 
to the lowest order in q±, the action (3) becomes Sll = So + Si, where 

So- paJ2(lPoM\xn\ 2 +yn6p*nXn) 

+ J dr(±5p* + ap e (d z x - 2 f^A z f + °f r%{dlxf 
Si * ^( 7 |O||s n| 2 +7'fi(s5n * 

Ho J dr (Wo) 2 + 2r^ 2 (s + B'^SA^ x z)f 



(5) 



(5)' 

where s T (the transverse component of s) is given 7 ' 8 ) by r B (d±x x z) with vortex spacing 
tb = ^/4>o/27vB, so the uniform displacement with vanishing q±, x a longitudinal phase 
variable, and Cqq the resulting shear modulus. Variational equations with replacement 
ifi — > u (real frequency) give eigenmodes of the vortex lattice. Note that the resulting 
dispersion ui = —i(a / ypo)~ 1 CQQ(q±rB) 4 for the shear mode 7 ' 8 ) at zero q z is different from 
that following from eq. (1), and hence that the extrapolation 1 ' 3 ) of the uniform (or equiv- 
alently, the single-) vortex dynamics to the interacting case with the shear mode is not 
justified within TDGL. Through the London limit (see the action (7) given below), the 
action (5) with the constraint s T = r B (d±x x z) is found to correctly give hydro dynamical 
results at longer distances than 



I = r B yC 66 r|/2apoCo> 

which is typically of the order tb- 

The absence of the dc uniform superfluid density (helicity modulus) is obvious from 
the expressions (4) and (5)', implying that a constant twist pitch 5A is transmuted into 
an uniform displacement of vortices keeping the free energy invariant. That is, as also 
indicated by Baym, 8 ) situation is different from an elastic matter 1 ) with free energy in- 
variant with respect to uniform displacements. Using (4) and (5)' and eliminating so, we 
easily obtain the mean field expressions on diagonal (i = j) and Hall (i ^ j) vortex flow 
conductivities 

M^ = 0, g2 )=p,o^ (1+r|e , 2)2 . (6) 
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Due to the constraint s T = r 2 3 (d±x x z ) imposed before introducing 6 A, the action (5), in 
which the degree of the positional order is directly reflected, does not give any nonlocal 
(and longitudinal) conductivity even when q± 7^0. It should be noticed that the nonlocality 
of the denominator in eq. (6), which results from the |<9 2 s| 2 term of eq. (5)' and is a typical 
example leading to a possible negative g 2 -term 5 ) in cr^ , is accompanied by a small length 
scale (~ t\b) insensitive to temperature and hence, is not relevant to the ordering of the 
system. Consistently with eq. (6), the nonlocal superfluid density, namely the real part of 
p S) ij(q, o; + i0 + ), becomes p a ^j {q±_ = 0, q z ) ~ p s o b~i 3 -r 2 B q 2 /2 in dc limit. As pointed out in 
ref.9, this q 2 behavior does not change through a 3D melting transition. 11 ) 

Next, let us here comment on eigen modes of the action (3) in 2D and nondisssipative 
(7 = 0) case, formally corresponding to the rotating superfluid 4 He at zero temperature 
described according to the Gross-Pitaevskii equation. Consistently with refs.12 and 13, 
the order parameter in eq. (3) will be divided into the amplitude p and phase (p; ift = 
v /pexp(i(/?) and, just as in the usual London approximation, the presence of the field- 
induced vortices will be taken into account through the topological condition 14 ) on d^(p = 
(d±<p, d T ip) by neglecting any fluctuation-induced vortices. As a result, we obtain the 
following harmonic action 

Sph = a IJ dT ^' 5pd ^ X + **o(0J-X " r- B \z x s)) 2 + ^-5p 2 + ^(disj^y (7) 
The only difference of this action from (5) (in 2D) is that the longitudinal current 

v s = 2ae (d ±X - r~ B 2 {z x s T )) (8) 

is nonzero in the action (7), which is necessary in considering longitudinal conductivity 
nonlocal in directions perpendicular to the field. When q±l > 1, the minimal coupling in 
eq. (8) between x an d s T is removed, and in principle the phase fluctuation behaves like in 
zero field case at such short lengths, while, if q±l < 1 and a gauge disturbance 5A±(q± ^ 0) 
is not substituted, v s ~ is established and thus, (7) reduces to (5). Variational equations 
resulting from the action (7) agree with those following from a dual representation used 
in ref.12, where a role of nonzero v s at nonzero frequency was stressed in a context of 
superfluid 4 He. The action (7) does not include the term 

S inc = i^| drY'^(sxd T s) z , (9) 

which is necessary in recovering the limit of incompressible fluid (i.e., dp = and v s = 0). 
This additional term with the coefficient 7", which should become 7', cannot be detected 
in this phase-only analysis for the action (3) (Hereafter, we will assume 7" = 7')- The 
eigenvalues following from .Sph+Sinc precisely coincide with those derived by Sonin 13 ) and 
give dispersions u 2 h ~ (cr^)" 2 bC 66 (q ± r B ) 4 / (1 + (q±l inc ) 2 ) and u 2 ~ (2^ 2 /7'r|) 2 (l + 
(q±knc) 2 ) of two modes corresponding to the shear (massless) and compression (massive) 
elastic modes, respectively, where Zi nc = r 2 3 ^/bpo/2a/$,o. The quadratic dispersion of cu s h 
at low q± is an origin of the destruction 7 ' 8 ) of the long-ranged phase coherence in the 



5 



vortex lattice. If Ti = 0, T 2 = l-nap^" B/4>0i an d v s = 0, the variational equation, with 
replacement —it — > t (real time), with respect to the shear displacement of <S p h + S- mc 
becomes the transverse part of eq. (1) with no external current. In this case the constraint 
v s = is different from that in the action (5) and rather corresponds to the limit of 
incompressible fluid. Following ref.13, this limit is valid only at shorter lengths than Z inc , 
which is at most of the order of several vortex spacings in a type II superconductor near 
the melting line in a moderate field. 

The above agreement with the well-known rotating superfluid hydrodynamics justifies 
the presence 7 ' 8 ^ of the minimal-coupling in v s between the phase field x an d shear dis- 
placement s T . Based on this, we assume below that the v 2 term appearing in the London 
limit (7) will be found at the static level by summing up many higher LLs in GL harmonic 
analysis around the mean field state. This should be expected as far as mean field results 
in London limit are recovered in GL theory. When combining this picture with eqs. (5) 
and (5)', we are naturally led to invoking the following action 15 ) appropriate to examining 
nonlocal conductivities: 

Q J 1 



PoaN(7lXgfi| 2 +%-f-|Sgd 2 ) + 



l s g nl 



(d ±X -r-*(zxs T ) - ^5A L ± f+r-^s L +B-\5Alxz)f (10) 
2 (po 

where 5A^ (<5A^) denotes the longitudinal (transverse) part of the gauge disturbance de- 
fined within the x-y plane. For later convenience, the shear modulus is assumed to have 
possible frequency and wavevector dependences, and the time scale of s was phenomeno- 
logically changed taking account of a possibility that it may have nonlocal corrections 
lq ~~ l(> 0)) due to an origin 1 ' 2 -' other than the freezing to the vortex lattice and irrele- 
vant to statics of vortex states at long distances. The amplitude mode leading to a 2 % 2 
term was neglected. The term (9) has to be included when considering a nonlocality of 
Hall conductivity. Again, the variational equation with respect to s T of the action (10) is 
different from that of eq. (1) with Ti = lnap^B / '<f> due to the presence of nonzero v s . 
Further, by neglecting the gauge disturbance S A and examining the dispersion of the shear 
mode in the vortex lattice, it is found that, due to the presence of the first term ~ |fi||x| 2 
in (10), setting v s = from the outset, as in eq. (1), always leads to an erroneous result 
on the dispersion. Rather, at low enough q, the constraint v s ~ is established as in eq. 
(5), and consistently, the second (dynamical) term |0||s T | 2 in (10) becomes unimportant 
compared to the first term. In the action (10), we assumed the coefficient of (d z s) 2 term 
to be the same as that in (5)'. Although this is not correct in the usual London limit, the 
details of this coefficient do not change main results in the next section. The first term in 
(10) having the same form as in the lowest LL case (5) is required within TDGL formalism 
and actually is justified, because the x variable to be related to the shear displacement 
at small but nonzero q± was shown 7 ) to, up to the lowest order in q± : become a phase 



+<oPo / / dr 
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change around ipo even if higher LLs are fully included. It is clear that this dynamical term 
associated with the Goldstone (shear) mode is found only by taking account of interactions 
among vortices. On the other hand, in the derivation in ref.3 of the eq.(l), the origin of 
contributions at nonzero frequency was attributed to variations of the order parameter in 
the vicinity of each vortex core which were assumed to be those of a single vortex motion. 
Consequently, the first (dynamical) term of the action (10) is overlooked in such an exten- 
sion of the single vortex dynamics to the vortex lattice. As is seen in the next section, due 
to the presence of this dynamical term, the destruction 7 ' 8 ) of the true off-diagonal long 
ranged order is reflected even in dc conductivities under a current perpendicular to B. 

§3. Nonlocal Conductivities 

It is straightforward, using the expressions (4) and (10), to find the nonlocal supercon- 
ducting contributions to dc conductivities. First, let us discuss the vortex lattice (solid) 
where C66(<? = 0, O = 0) ^ 0. In this case the longitudinal (|| qj_) part of the diagonal 
conductivity a xx and the Hall conductivity are given by 

is)L(a\-n r B %Qi + lQlr 2 B [q 2 z + m±/rB) 2 ] 2 /2 
[Q) Ps0 2e fe 2 2 (l + S 2 r|/2)+;Vsi] 2 ' 

2 / 2 

a ^ {q) = Ps0 2# (l + q ^ B /2M(l+ Z q ^rlm + Q 2 Q 2 J 2 Y 

using the length / oc y/C^ defined in § 2. In a layered material with mass anisotropy 
M/m(> 1) and under a field perpendicular to the layers, q 2 appearing in expressions (11) 
is replaced by 2(1 — cos(q z s))m/ M s 2 , where s is the layer spacing. The transverse part a xx 
of the diagonal conductivity merely becomes, even in the liquid regime, 7 g a xx (q± = 0, q z )/~f 
(using eq.(6)) , and hence, in contrast to (11), is not affected by the positional correlation 16 ) 
and not relevant to the channel flow situation, proposed in ref.l, where rather a xx is 
measured. We note that the combination ~ q 2 + qj_l 2 seen in denominators of (11) is the 
dispersion of the 3D shear (Goldstone) mode destroying the true off-diagonal long ranged 
order. 7 ' 8 ) In general, the magnitudes of the conductivities (11) significantly depend on the 

relative size of \q z \ and \q±\. For instance, when \q z \ > l\q±\/rs ~ both a x s x L and a x s J 
are well approximated by cr^- ( expressions (6)) with replacement 7 — > 7 g . On the contrary, 

when \q z \ < lq\, a x s x L typically becomes the value at zero q z (or in 2D case); 7/9 s o/ '<Zi_£o> 
which is the same as the 2D result in zero field. 17 ) Although this result has also been 
suggested in ref.5, it is unclear to us whether the authors in ref.5 have distinguished the 
case \q z \ > \q± \ from the case \q z \ < \q±\. This divergent behavior is a reflection of the fact 
that, in this case, the conductivity is primarily determined by the phase variation x an d 
not by the so-called vortex motions d t s = £? _1 (E x z) implicit in the last term of the action 
(5)', where E is the external electric field. Correspondingly, a x s J decreases like ~ (qz/q±) 2 
and vanishes, for arbitrary (nonzero) q±, at zero q z or in 2D case. In other words, as 
the nonlocality parallel to B becomes negligible comparted to that perpendicular to B, 
contributions of a finite ^ q are covered and overcome by the perfect positional correlation 
(i.e., the infinite shear viscosity) making the Goldstone mode well-defined. Particularly, 
these features in \q±_\ > \q z \ are quite remarkable, because, as mentioned in relation to 



7 



the actions (5) and (5)', their origin (i.e., a nonvanishing v s ) is of higher order in the 
wavevector q± and hence, is usually neglected. On the other hand, these features are lost, 
at fixed q± and q z , as / becomes shorter, i.e., with increasing field. It is consistent with the 
absence of the nonlocal conductivity in the vortex lattice constructed within the lowest LL 
(see a sentence following eq. (6)). Unfortunately, we cannot predict precisely a threshold 
field below which these nonlocal effects are measurable. In addition, note that, as seen 
above in q z — > limit, the infinite shear viscosity of the vortex lattice does not mean a zero 
value of the nonlocal resistivity. We emphasize that the remarkable difference between the 
uniform result (6) and this large conductivity in the case with zero q z and small q±(^ 0) 
is a consequence of the fact that, at zero q±, the shear displacement s T decouples with the 
zero mode x an d changes into the uniform displacement. As is seen by comparing (10) 
with (5)', the ac conductivities in the vortex lattice, when q± — > 0, reduces to the vortex 
flow results (6) except corrections irrelevant to positional correlations. This is a typical 
example of a crucial difference between ac and dc responses in a kind of ordered state. 

(§)£, _o ( s ) 

The above results that a xx ~ q ± and a xy = when q z = and q± f 1 are valid 
even for a state with nonzero C§§{q±_ ^ 0, |0| = 0) such as the 2D hexatic liquid state 
which may be possible in low enough fields (The 3D hexatic phase invoked in refs.l and 
10 is inconsistent with the first order freezing to the vortex lattice and never realized 18 -* 
in real superconductors, and hence need not to be considered). They are also applicable 
to real systems with pinnings and hence with a finite correlation length £/,, of the bond 
orientational order, if q±^h > 1- Therefore, it is valuable to examine these dc linear 
responses in real experiments in relation to determinations of the position of the freezing 
to the vortex lattice and of the existence or absence of the hexatic phase in 2D systems. 

(§)£, o 

In particular, it should be noted that the behavior a xx ~ q ± at nonzero q± as well as 

the vanishing a x s y implies p xy ~ P xx °~x y ,Ni where p xx (p xy ) is the total diagonal (Hall) 
nonlocal resistivity, and <J xy ,N is the extrapolated normal Hall conductivity. 

When trying to understand the liquid regime (disordered phase) in the present ap- 
proach, some comments are necessary. As shown in ref.9, the harmonic analysis cannot 
be used in the situation where the contributions with nonzero reciprocal lattice vectors 
are negligible even if the amplitude fluctuation of ip around i/jq is negligible. Even in such 
situations the nonvanishing transverse diamagnetic susceptibility results from the static 
superconducting fluctuation 6 ) and, in layered systems, shows a dimensional crossover 9 ) 
due to a competition between the layer spacing and a finite phase coherence length paral- 
lel to B which cannot be detected in the harmonic analysis. Consistently, this dimensional 
crossover will be seen also in the nonlocal conductivity. Unfortunately, it is practically 
difficult at present to find possible nonlocal corrections, associated with the freezing to 
the vortex lattice, to the vortex flow expression (6) according to the fluctuation theory, 6 ) 
but the mean field (i.e., harmonic) analysis may become a guideline to an extension of 
the analysis in ref.6. In addition, it is possible that the present analysis based on the 
phase-only model (i.e., based on the mean field solution) will provide an essentially correct 
result at relatively short scales in directions perpendicular to B in the disordered state. 9 ) 
For these reasons, we will focus on the region just above the freezing point to the vortex 
lattice and only consider the nonlocal conductivities with vanishing q z (see, however, § 
4), in which case the dimensional crossover does not appear, and hence, the use of the 
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harmonic approximation may be permitted. 

In this case, the corresponding results to (11) are given by 



„{i)L (n .„ r B % + v^b/po^ 

XX ^^-^-P^2H l + 5ir|(7 ff + ^r|/poa)/27' 

2 / 

" S(4X:4 ' = 0) = ""it 1 + Sft + 7 W ir|/p a)/2r (12) 

where the Maxwell form 19 ) C66(<?, ~ 77 1 1 with shear viscosity 77 was assumed 20 ) 
for the dynamical shear modulus. Interestingly, the expressions (12) smoothly lead to 
the results (11) for the vortex lattice at zero q z by taking 77 to be infinity. The terms 
1 + ^(qxjB) /^Po® m the denominators of expressions (12) again originates, through 
a nonzero v s (i.e., a higher spatial gradient), from the spectrum of the shear mode 
u> ~ — iC , 66(<?±^s) 4 , which cannot be detected in eq. (1), while the nonlocality in the 
numerator of o~ x l l L , resulting from nonlocalities of the time scale for the displacement field, 
is the same as that expected from (1). As a result, at large q± of the order r^ 1 , both (11) 
and (12) are well approximated by the usual one (6), suggesting that, for such a rapid 
variation of the external current, not only collective effects but also the sample disorder 
(pinning) is irrelevant as far as the time scale to be affected by the sample disorder is % 
and not 7. In general, in the case with a pinning effect enhancing 7^, the viscous effect ac- 
companying the freezing to the vortex lattice becomes negligible. 1 -* Further, the expressions 
(12) suggest that the low temperature limits of the nonlocal conductivities with nonzero 
q±_ in 2D systems with pinnings, where we have no transitions at nonzero temperatures, 
again become the corresponding, and above-mentioned, results (independent of <y q ) in the 
vortex lattice, again leading to the relation p xy ~ P xx o~ xy ,N- This relation should also 
be observed in a clean 3D-like situation with vanishing q z and a large 77. An apparently 
similar relation was found 21 ) in the pinning-dominated region of BSCCO. 

§4. Discussions and Conclusion 

We will discuss consequences of possible viscous effects 1 ' 2 -* deep in the liquid regime 
of pinning free systems, other than 77, which should appear through nonlocalities of 7 5 . 
According to eqs. (11) and (12), when \q z \ <C \q±_\, it is difficult to practically divide, in 
Cx2 L , contributions of 7 g from the positional correlation effects, and hence, we will only 
consider the case with vanishing q±_ but nonzero q z , where the diagonal conductivities can 
be always expressed by (6) with replacement 7 — > 7 g . Further we will focus on the 3D 
region below the dimensional crossover 9 -* (We note that this dimensional crossover, formally 
equivalent to that resulting from the entanglement picture, 1 ) has nothing to do with that 
argued through experiments in ref.4). For this case, it is often argued 22 ) that thermally 
activated cutting (and reconnection) processes among vortices induced by a nonuniform 
(Qz 7^ 0) current will lead to a very long length scale and significantly increase 7 g even if 
the (thermally-induced) entanglement is absent. 2 ) This picture seems to suggest that the 
time scale grows unlimitedly even in the vortex lattice upon cooling, and hence that the 
vortices bent by the current cannot move at low enough temperatures. It should be noted 
that such an argument would also be applicable to a nondissipative case (with zero 7 but 
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nonzero 7') by imagining a nonlocal growth of 7'. However, it is even unclear to us if 
this is a correct argument in the context of linear hydrodynamics. In superfluid 4 He, the 
reconnection process does not 23 ) seem to need a remarkable slow dynamics. In addition, 
the algebraic system-size dependence of an onset temperature 4 ' 24 -* of apparently nonlocal 
vortex motions is inconsistent with the argument based on the thermally activated vortex 
cutting processes in which the temperature variation is dominated by the Arrenius factor, 
because the measured algebraic size dependence inevitably means an algebraic decrease 
of the cutting barrier itself with increasing the system-size, although intuitively such a 
decrease of the barrier should not be expected. Further, the resistivity data 24 ) parallel 
to B in a heavily twinned sample seem to intuitively contradict those in a more 3D-like 
situation. 25 ) Clearly, experiments in twin-free samples are necessary. 26 ) Theoretically, there 
seems to be no acceptable reasons why such a nonlocal growth of 7 g has to be expected. 
Firstly, our analysis shows that, in contrast to the phenomenology in ref.2 and 5, viscous 
terms possible in 7^ cannot reduce to elastic terms in the vortex lattice. Secondly, the (if 
any) entanglement 1 making a growth of 7 9 possible must disappear at least deep in the 
vortex lattice. 27 ) In ref.9, the presence of a large and positive nonlocal (~ g 2 ) contribution 
in 7j was questioned on the basis of the observation that, in contrast to the pinning- 
induced activation form for the time scale, any nonlocal (i.e., (/^-dependent) growth of the 
time scale may be incompatible with the uniform vortex flow due to the nonlinearity (i.e., 
summations with respect to internal <7 z 's) in Kubo formula in the nonGaussian fluctuation 
theory. 6 ) At the present stage, however, this opinion is not still conclusive. 

As already commented in the context, the present analysis cannot clarify a charac- 
teristic field below which the positional correlation-induced nonlocal responses studied in 
this paper and the resulting relation p xy ~ p xx become measurable. In order to estimate 
this theoretically, it is necessary to know how a nonzero v s , required through the agree- 
ment with the rotating superfluid hydrodynamics, is found consistently with the action (9). 
From the fluctuation theory, 6 ) the nonlocality in the numerator of a xx L may be expected 
to be detected as a consequence of renormalizations of the time scale for vortex motions. 
On the other hand, it is extremely difficult to detect nonlocalities in the denominators of 
conductivities (12), because their origin consists in the details of physics at short scales 
which, as mentioned in ref.9, cannot easily be captured in the fluctuation theory from 
higher temperatures. 

In conclusion, a hydrodynamics for vortex states in type II superconductors consis- 
tent with the rotating superfluid hydrodynamics has been presented in order to study the 
nonlocal conductivities perpendicular to the magnetic field. In particular, the conductiv- 
ities with nonlocality only in perpendicular directions to the field and in a situation with 
large shear viscosity are essentially the same as those in zero field case, and the vanishing 
superconducting part of the Hall conductivity under such a situation will be useful in ex- 
perimentally judging the existence or absence of the 2D hexatic phase in type II supercon- 
ductors. It is interesting to experimentally examine the predicted relation p xy ~ P xx &xy,N 
not only in dirty samples but in a clean sample such as untwinned YBCO. 
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